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We perform a self-consistent relativistic QRPA (RQRPA) calculation of 2vj3j3-decay based on 
relativistic BCS (RBCS) mean field theory results for odd-odd intermediate nuclei ^^Sc, ^^As, 
82gj.^ '"''Tc, '^^I, and ^^''l. The RBCS equations that resemble the non-relativistic ones are con- 
structed from a Dirac-Gorkov variational functional. We use the parameter set NL\ for the (7, (O 
and p mesons. The RQRPA equations are solved for the residual n + p interaction by employing 
the same parameters used in the RBCS for the latter meson, and experimental values for the pion 
and nucleon. The RQRPA results for the 2v/3/3 matrix elements are similar to those obtained 
within the QRPA and the shell model. 
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1. Introduction 

In nature there are about 50 nuclear systems in which the single j8 decay is energetically 
forbidden, and jSjS decay turns out to be the only possible mode of disintegration. It is the nuclear 
pairing force which causes such an "anomaly", by making the mass of the odd-odd isobar, (A^ — 
1,Z + 1), greater than the masses of its even-even neighbors, {N,Z) and {N — 2,Z — 2). The modes 
by which this decay can take place are connected with the neutrino (v)-antineutrino (v) distinction. 
In fact, they are defined by the transitions: 

p + e^ + Vrh, (1.1) 

the neutrino V being left-handed (LH) and the antineutrino v right-handed (RH) because of parity 
non-conservation in weak interactions. Therefore, regardless of the Dirac (v / v) or Majorana 
(v = v) nature of the neutrino and independently of conservation of helicity, the two-neutrino 
mode (2vj8j8) decay can occur by two successive single jS-decays: 

{N,Z) ^ (A^-l,Z+l) + e- + v 

(A^-2,Z + 2) + 2e" + 2v (1.2) 

passing through the intermediate virtual states of the (A'^ — 1 , Z + 1 ) nucleus. Yet, the occurrence of 
the neutrinoless jSjS decay (OvjS/3): 

(A^,Z) ^ (A^-2,Z + 2) + 2e~ (1.3) 



is much more convoluted since the right-handed neutrino emitted in the first step of (1.2) has the 



wrong helicity to be reabsorbed in a second step and to give rise to (1.3). For massless neutrinos 
there is no mixture of left and right handedness, and the OvjSjS -decay cannot occur, regardless 
of the Dirac or Majorana nature of the neutrino. Yet, experiments with solar, atmospheric and 
reactor neutrinos have provided remarkable evidence in recent years for the existence of neutrino 
oscillations driven by nonzero neutrino masses and neutrino mixing [|], ^]. Once the neutrino 
becomes massive, the helicity is no longer a good quantum number. Then, if the neutrino is in 
addition a Majorana particle with an effective mass {my), the mixture of Vlh in Vrh is proportional 
to {my) /Ey, and OVj8j8 -decay is allowed. ^ This fact inspired experimental searches in many nuclei, 
not only for the Ovj3j3 -decay but also for the 2vj3j8-decay, since these two modes of disintegration 
are related through nuclear structure effects. In fact, their half-lives can be cast in the form: 

T^y ='^lVy^2y, T^J =%y^^y{my)^, (1.4) 

where '^'s are geometrical phase space factors, and the ^'s are nuclear matrix elements (NME's). 

and present many similar features, to the extent that it is frequently stated that we shaU 
not understand Ovj3j3 -decay until we understand 2vj8j8 -decay. 

Quantum hadrodynamics (QHD) aims to describe the nuclear many-body system in terms of 
nucleons and mesons |^. Proposed initially as a fully renormalizable quantum field theory, at 

'For simplicity we assume that right-hianded weak currents do not play an essential role in the neutrinoless decay. 
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present it is seen as an effective field theory, derivable, in principle, from quantum chromodynam- 
ics [Q]. Relativistic mean field theory (RMFT), which can be thought as a mean field (Hartree) 
approximation to QHD, has been applied with great success during the last decades to account for 
nuclear matter saturation and the ground state properties of finite nuclei along the whole periodic 
table. Through a relativistic version of the random phase approximation (RRPA), various excited 
states and the Gamow-Teller (GT) and Fermi (F) resonances have been studied in this context as 
well [^, 1^ ^, |8|, ^, 1^, 11, 12, 13]. When this approximation is based on the Hartree-Bogoliubov 
(HB) or BCS approximation, it is called the relativistic quasiparticle RPA (RQRPA). This approach 
has been used to evaluate several weak interaction processes, such as beta-decays, neutrino-nucleus 



reactions, and muon captures [ 14, [T^, 16]. Here the first application of the RQRPA to 2vj8j3 -decay 
is made. 



2. 2v/3/3 matrix element 

Independently of the nuclear model used and when only allowed transitions are considered, 
the 2vj8j8 matrix element for the |0t) final state reads [17] 



^2v(/) = I (-)'I 



A=0.1 



(o;ii<iiA+)wii<iio+) 



-^■^!v{f) + -^2v\f) (2-1) 



where Eq and Eq+ are, respectively, the energy of the initial state |0+) and of the final states |0l 
The summation goes over all intermediate virtual states and 



(2A + l)-""i;{p||Oi||n) (4c,,),, with 



pn 



00 = 1 for F 

01 = a forGT 



(2.2) 



are operators for the j3 -decay. The corresponding j8+ -decay operators are ~ i^x ^ ■ ^^'^ 



total j8* strengths 



5f = (2A + l)-i£|(A+||<||0+)p, (2.3) 



obey the well-known single-charge-exchange Ikeda sum rule (ISR) ^ for both the F and the GT 
transitions: 



-i/n+ir/?/5' 



'X '^A 



Similarly, the jSjS -decay strengths 



(2A + 1)-^£K0 



"]o|0+)=A^-Z. 



|0+)P 



(2.4) 



(2.5) 



obey the double-charge-exchange sum rules (DSR): 

5f =5f"-5f" = (2A + l)-i(0 



(2.6) 
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which when evaluated give [19|: 



5^ 



2{N-Z){N-Z-\), 
2{N-Z) (n-Z-\+2S^^ 



(2.7) 



where C is a relatively small quantity and is given by [ ]19| , Eq. (5)]. The DSR are as important to 
/3j8 -decay as the ISR is for simple jS -decay. 

Contributions from the first-forbidden operators appearing in the multipole expansion of the 
weak Hamiltonian, as well as those from the weak-magnetism term and other second order correc- 
tions on the allowed 2vj3j3-decay, have been examined rather thoroughly [20, 



3. Charge-exchange QRPA 



The pn-QRPA was formulated and applied to the allowed jS ^-decays and to the collective GT 



resonance by Halbleib and Sorensen (HS) in 1967 [22]. They solved the QRPA equation 



A B 
B A 



X 



(On 



X 



within the pn quasiparticle (qp) space for the BCS vacuum 

|0^) = Ylit^P+Vpclcl)Y[iUn + Vnclcl)\), 



(3.1) 



(3.2) 



of the initial nucleus {N,Z), where |) stands for the particle vacuum. The transition matrix elements 



are: 









= 1 


^pVn'^pn;l 








pn 








l|0+) 


= 1 


VpUnXpn-l 



(3.3) 



pn 



the ISR ( p3[ ) yields N — Z, and the ground state correlations (GSC) in ( p^D play an essential role 
in suppressing j8+ -decay. 

Intensive applications of the QRPA to jSjS -decay began only about 20 years later when Vogel 
and Zirnbauer [ p3| ] discovered that the j3+-decay suppression mechanism could also be invoked 
to explain the quenching of the 2vj8j3 decay rates. Their adaptation of the HS model in essence 
implies: 1) To use a second BCS vacuum for the final nucleus {N — 2,Z + 2), and to solve a second 
QRPA equation for the intermediate 1^ states with the ISR equal toN — Z — 4, and 2) To substitute 
(|l|)by theansatz 

^2v 



21 

aa' 



(i«'ii<iio^)(C^iu)(uii^riio+) 

(0,+ 



(3.4) 



To circumvent the nonphysical averaging procedure implicit in the overlap (l^/lla), a different 
recipe for the application of the QRPA to the jSjS -decay has been introduced [25, 26, 27, 28, 29, 3^], 
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which continues to involve two BCS ground states but deals with only one set of QRPA solutions. 
Finally, the procedure has been simplified even more by performing a straightforward adaptation 
to the jSjS-decay of Cha's prescription for the evaluation of single j3-decay [ ^T] ] within the HS 
model, which implies solving both the BCS and QRPA equations for the intermediate {N —\,Z-\-\) 
nucleus. The ISR then gives N — Z — 2, and the above expression becomes [32] 



(\+\ 



"||0+)(1 



110+) 



(3.5) 



Numerical tests show that all the approaches above furnish quite similar results for 1 25, 26, 



28, 29, 32]. Thus, for the sake of simplicity, we make use of the last expression in the present work. 
When applied to the j3j3 -decay, the QRPA turns out to be an uncomplete model, since it deals 



with and 2 qp states only, while to evaluate Eq. (1.2) it is necessary to consider at least up to 4 



qp states. Moreover, the QRPA can say nothing regarding the DSR given by (2.7), nor it can be 
used to describe the )8 )8 -decays to excited final states. How this can be implement is explained in 
Ref. (IT; 



4. Relativistic charge-exchange RQRPA 

The RQRPA is based on the relativistic HE (RHB) approximation for the RMFT. It was formu- 
lated in Ref. [34] for charge-conserving excitations, and extended to charge-exchange excitations 
in Ref. [35]. In the present work we approximate the RHB equations and put them in a form that 
resembles the non-relativistic BCS equations. To do this we start from the variational functional 



{(Or + ^r)d{x-y)-h,{x,y) Al{x,y) \ ( U^ix) 

A,(x,y) {(Ot-Ht)5{x-y)-ht{x,y) \roV,{x) 



(4.1) 



based on the Dirac-Gorkov equation [ p3| , Eq.(39)] with notation t = p or n. Here Ut{x) and Vf(x) 
are the normal and time-reversed Dirac spinors corresponding to solutions of this equation with 
positive and negative-frequency ft)^. The Lagrange multipliers /if are determined by requiring that 
the expectation values of the baryon number operators yield the desired values of Z and A'^. Dirac 



Hamiltonian operators ht{x,y) and pairing fields At{x,y) are given by [33, Eqs. (40) and (49)]. 
Next we use the anzatz 



U,ix) \ _^ / ut%{x) 
YoVtix) ) i Vf^,(x) 



(4.2) 



where and Vt sue numbers, and '^t (x) are the Hartree mean-field wave functions, satisfying the 
equation / d^yht{x,y)^/^t{y) = £r^r(x), with Et being the single-particle energies. After performing 
this replacement one maximizes ( |4. 1[ ) with respect to the coefficients Uf, obtaining in this way the 
relativistic BCS (RBCS) equations for m, and v,, similar to the non-relativistic ones. 

Since the pion does not participate in the RBCS, the Lagrangian density is determined once 
the masses of the nucleon and mesons a, CO and p, the coupling constants of mesons with the 
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nucleon, and the self-interaction constants of the meson a, g2 and g3 are given. Several sets for 
these parameters are known in the literature. Here we use the NL\ set[^]. 

We solve the RBCS and the Klein-Gordon equations numerically by expanding the mesons 
fields and the fermions wave functions in complete sets of eigenfunctions of harmonic oscillator 
(HO) potentials. In actual calculations, the expansion is truncated at a finite number of major shells, 
with the quantum number of the last included shell denoted by A^^ {Nb) for fermions (bosons). The 
maximum values are selected so as to assure the physical significance of the results. The oscillator 
frequency for fermions is given by HcOq = AlA^^I^ MeV and the maximum number of oscillator 
shells for fermions and bosons is given by Nf =Nb = 20. The Coulomb field is calculated directly 



in configuration space. The same procedure was used by Ghambir et al [37] in their approach to 
the relativistic mean field. 



The RQRPA equations ( |3.1[ ) are solved by employing for the residual interaction the same 
parameters used in the RMFT to obtain the discrete basis of qp states within the RBCS approx- 
imation. Yet, in dealing with isovector excitations it is essential to include, together with the p 
meson, the n meson as well ^. Here, the experimental values of the pseudoscalar pion-nucleon 
coupling and the pion mass were used, i.e., = 1.00, and m-n = 138.0 MeV [^, |6|]. Since Fock 
terms are ignored in the RMFT, for the sake of self-consistency we must omit the exchange matrix 
element of the residual interaction V = + Vp in the sub-matrices A and B. 



5. Results 



Table 1: Matrix elements in units of (MeV) The experimental values are obtained from the half- 
lives T2v compiled by Barabash pq] and the ^#2v values from Ref. [B9| for the bare axial-vector coupling 



constant ; 



= 1.25. 



j3j3 -Decay 



EXP SM® SMren[||] QRPA [||] RQRPA 



>48Ti 

^^^Se 
^82Kr 



48Ca ^48sc 
ioomo^iooTc^iooRu 

128xe^l28l^l28xe 

I30xe^i30i^i30xe 



0.049 ±0.003 
0.140 ±0.005 
0.098 ±0.004 
0.239 ±0.007 
0.049 ±0.006 
0.034 ±0.003 



0.047 
0.116 
0.126 

0.059 
0.043 



0.055 
0.206 
0.224 

0.116 
0.085 



0.058 
0.064 
0.077 
0.065 
0.076 
0.061 



0.066 
0.055 
0.088 
0.062 
0.076 
0.070 



In Table |l| the experimental values of the matrix elements are compared with values 
obtained from several calculations. The later critically depend on the adopted value for the effective 



axial- vector coupling constant ^a- The QRPA results [[29|], as well as the present ones correspond 
to = 1> i-^-, to a quenching factor of ^ = ^a/^a = 0-8 when the bare value is = 1.25 [^]. 
Bearing in mind that only a very tiny fraction (< 0.1%) of the sum rule 5^^ goes into the final 
state = |[T^], we can say that both agree reasonably well with experiment. One should 
also remembered that the QRPA calculations have been monitored by the restoration of the SU (4) 
symmetry while the RQRPA were not. In the shell model (SM) study pO| ] different q values were 
used in different nuclei, namely q = 0.74 in "^^Ca, q = 0.60 in ^^Ge, and ^^Se, and q = 0.57 in ^^^Te, 
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and ^^°Te. These results are listed in the third column of Table |TJ For the sake of comparison, the 
same results renormalized to ^ = 0.8 are shown in the fourth column (labelled as "SMren"). 

6. Conclusion 

A variational functional based on the Dirac-Gorkov equation is used to obtain the RHB equa- 



tions in the form of the non-relativistic BCS equations. The RQRPA equations (3.1) are solved 
for the residual 71 + p interaction by employing for the latter meson the same parameters used in 
the RMFT The RQRRA^ results for the 2vj3j3 matrix elements are of the same order magnitude 
as those obtained within the QRPA and the SM. Bearing in mind the small fraction of double GT 
strength going to the 0+ final state compared with the GT DSR, as well as the uncertainty involved 
in the quenching factor q, it is difficult to discern which of the three calculations is better and 
which is worse. Despite this we are planning to apply the RQRPA model to study the neutrinoless 



j8j8 -decays using the formalism developed in Ref. |21] as well. 
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